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We numerically solve the non-linear Poisson-Boltzmann equation for two cylinders confined by 
two parallel charged plates. The repulsive electrical double layer component of the cylinder pair 
potential is substantially reduced by confinement between like-charged plates. While the effective 
cylinder surface charge is increased by the confinement, the effective interaction screening length 
is reduced, this effect being dominant so that the repulsive confined cylinder- cylinder interaction 
potential is reduced. 



o 

03 



T3 

o 
o 



> 

(N 
O 
(N 
ON 

o 

OO 

S 

i 

C 

O 

o 



X 



I. INTRODUCTION 

Recent experiments have cast doubts that the DLVO 
pair potential is correctly describing the pair interaction 
between like-charged colloids in aqueous suspension in 
a confined region where the colloid motions are being 
restricted by the confining double layers. Long range 
attractive potentials of order 1 KbT in strength have 
been observed Additionally there is interest in 

whether shorter range interactions between like-charged 
cylinders in monovalent electrolyte can become attractive 
under certain circumstances |]-|^] ■ Here we examine the 
problem of long range interactions between parallel like- 
charged cylinders confined between like-charged plates by 
numerically solving the two dimensional, confined, non- 
linear Poisson-Boltzmann equation. 

II. NUMERICALLY SOLVING THE 2D, 
CONFINED, NON-LINEAR 
POISSON-BOLTZMANN 
EQUATION 

Two water-solvated like-charged cylinders experience 
an electrical double layer (EDL) repulsion, their behav- 
ior governed by the non-linear Poisson-Boltzmann (PB) 
equation, which dictates both the potential and simple 
ion concentration distributions in their vicinity: 



k 2 sinh i 



(1) 



Where V 2 is the scalar Laplacian operating upon the 
dimensionless potential <f>, which equals the scalar poten- 
tial tp divided by KeT/ze, with e being the quantum of 
charge, and z the charge of a single counterion. k _1 is 
the Debye screening length. Assuming the simple ions to 
be monovalent K B T/ze= 25.69mV at 298 K. 

The EDL interaction forms the repulsive component 
of the DLVO potential between two like-charged col- 
loidal particles. At constant thermodynamic volume the 
Helmholtz free energy is appropriate for describing the 
EDL interaction between two like-charged colloidal par- 
ticles. At constant surface potential, the EDL interaction 
energy has three parts: an attractive electrical term (neg- 
ative cylinder surface residues attracted to the positive 



charge counterion cloud between them), a repulsive en- 
tropy term (ion/solvent configurational entropy decreas- 
ing with decreasing cylinder separation), and a repulsive 
chemical potential term (counterion number is decreased 
by surface charge condensation as the cylinders approach 
each other). While at constant cylinder surface charge, 
the EDL has only two parts: an attractive electrical term 
and a dominant repulsive entropy term (with constant 
surface charge, counterion number is constant). 

The two dimensional problem of circles confined by 
line charges is equivalent to parallel cylinders confined 
by walls in three dimensions. Numerical computations 
of the nonlinear Poisson-Boltzmann (PB) equation were 
performed principally for two circles of constant (dimen- 
sionless) charge density (a — d ^ d %" where n is the normal 
to the surface) confined by two constant (dimensionless) 
potential line charges (4>n nc ) as sketched in Fig. ||. The 
boundary condition was usually constant charge on the 
cylinders and always constant potential on the confining 
charged walls because we supposed the cylinders to pos- 
sess strong acid groups, typical of surface groups such 
as polystyrene sulfonate, while we assumed the confining 
walls were made of glass, which contains a high density 
of weak acid silanol groups. 

930 





t ! 




1 c 



1 



FIG. 1. Geometry for solving the 2d, confined, non-linear 
Poisson-Boltzmann equation. Lengths are in grids of a 2d 
square lattice, where the reservoir, or bulk screening length is 
6-30 grids. Typically, a constant potential boundary condition 
(b.c.) is imposed on the confining line charges, a constant 
charge (potential slope) b.c. is imposed on the circles, and a 
free b.c. is employed at the left and right boundaries of the 
space (extrapolating the local potential slope, which equals 
zero far from the circles). The space's width was chosen so 
that the circle-circle interaction potentials were insensitive to 
increases in width, c is the closest distance between the circle 
and the line, and r is the center-center separation between 
the circles, all in grids. 

Our use of a constant potential glass boundary con- 
dition needs further clarification. This is a standard low 
surface potential, weak acid surface group boundary con- 
dition @ and is enforced at the outer Helmholz plane 
(OHP) where the compact, or Stern layer ends, and by 
association dissociation equilibria of weak acid surface 
groups and also two dimensional mobility of counterions 
within the compact layer. This layer is a highly concen- 
trated monolayer of aqueous counterions typically con- 
taining 90 percent of the glass countercharge and well 
over 100 millivolts of potential change so that the result- 
ing glass boundary condition which faces the electrolyte 
at the OHP is a low (usually well less than 100 millivolt) 
constant potential pf]| , p7t . Our results depend upon this 
glass regulating constant potential boundary condition's 
ability to stand up under compression. Should it fail then 
there would follow electric field lines leaking into the low 
dielectric glass and polarizing electron clouds. It would 
then become important to account for unscreened inter- 
actions between these polarized electron clouds. This 
paper discusses the limited constant potential confining 
boundary case. A complete treatment would include the 
field inside the glass, model the glass- water interface with 
a regularizing boundary condition flq| , as well as calcu- 
lating the field between the charged circles. Although it 
is expected that the constant potential glass boundary 
condition would remain valid at compressions between 
surfaces down to one screening length |l5| there is some 
evidence that under extreme compression it is in fact the 
case that field lines are traversing through the external 
dielectri c |3j ] . The inclusion of regularizing boundary con- 
ditions [ |16| is a direction for future research. 

There is zero electric field within the circles. Note 
that the constant surface charge boundary condition on 
the circles will not be a problem because field line pene- 
tration into the low dielectric circles causes only a small 
perturbation in the circle-circle interaction (see Fig. 9). 

A periodic boundary condition was employed at the 
left and right edges of the space (at these edges we extrap- 
olated the local potential slope in the direction parallel to 
the confining lines, which was equal to zero, since these 
boundaries were many tens of screening lengths separated 
from the circle edges; see Fig. 4). The space's width was 
chosen so that the circle-circle interaction potentials were 



insensitive to increases in this width. 

The geometrical details of our numerical simulation are 
as follows. Our circle's radius a was fixed at 36 grids of a 
2d square lattice, the screening length kT 1 varied from 6 
to 30 grids, while the space was made 930 grids in length 
in order to avoid edge effects. The distance from the cir- 
cle's surface to the line charge c was varied from 6 to 
\ screening lengths. Venturing much below kc = \ was 
found to introduce coarse graining errors, and a more so- 
phisticated multi-grid relaxation algorithm jl2| would be 
necessary in order to accurately track the highly curved 
potential function in this very strongly confined regime. 
Potential fields for fixed circle and line boundary condi- 
tions, c, and k _1 were obtained by numerically solving 
the non-linear Poisson-Boltzmann equation on a Silicon 
Graphics work station having a MIPS R4400 processor, 
and by using a successive overrelaxation (SOR) |lj] al- 
gorithm, employing an SOR factor of 1.85. Total free 
energies TFE for a given separation r were obtained via 
the Helmholtz prescription from Overbeek so that 
a circle-circle interaction potential could subsequently be 
constructed. Then we made a two parameter fit of the 
Helmholtz interaction energy to 



V EDL3B =*GF<L 



(2) 



4 being the cylinder's charge per unit length. The fit 
netted K~j~o C ai an d a geometric prefactor GF depending 
upon local screening length and particle radius. Between 
spheres in three dimensions |l3| 

„K,i oca i2a 

GF 3D = —— rj, (3) 

(1 + KlocaiaY 

while between cylinders in two dimensions 

K iocai 2a 

GF 2D = — — ^. (4) 

(1 + 2K loca ia) z 

These fits revealed that use of the bulk screening length 
was not appropriate for the confined circle-circle in- 
teraction. As the circles were increasingly confined by 
the line charges the circle-circle interaction's effective 
screening length was found to decrease in a systematic 
way, with the line charge counterions screening the circle- 
circle interaction. In addition, the circle's effective charge 

VGF was found to be increased by the confinement. 

Overbeek's Free Energy bookkeeping |ll[] goes as fol- 
lows: 

Total Free Energy = TFE = EFE + CFE (5) 
EFE = Electrical Free Energy = EE - TAS (6) 
CFE* — Chemical Free Energy* = 

= -2 / ((V0) 2 + K 2 (/.sinh0)dA (7) 

J A 
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EE* = Electrical Energy* = / (V0) 2 eL4 (8) 

J A 

TAS* = -2k 2 / (1 - cosh 4> + 4> sinh 4>)dA (9) 
J A 

Where * indicates a dimensionless 2d energy, dS repre- 
sents a small element of constant potential boundary, and 
dA a small area of electrolyte. Converting a 2d dimen- 
sionless free energy into a 3d dimcnsionfull one requires 
multiplying by keT, by and by one half the num- 
ber of Bjerrum lengths the circular cylinders extend into 
the z direction 113]. Take Fig. || as an example: if its 
cylinders were 500nm in radius, K~~ =280nm (in order 
that Ka— 1.8), and the cylinders were lOOOnm in length, 
then we should multiply their dimensionless 2d energies 
by k B T-J-±220iiHi = 57 k T i e . at Kr= io, KC =6 the 

J D 47T 1.4nm ° 1 ' 

cylinders would experience about a 0.6 keT repulsion. 
The free energies obtained from numerical solutions of 
the PB for interacting flat plates compared well against 
the tabulated values in the Verwey-Overbeek monograph 
p7[ , and also against Israelachvili's figure 12.10 |lE{ ]. 
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FIG. 2. Reduction of interparticle repulsion between a 
pair of charged rods due to the rods being confined by a pair 
of charged plates. 2d dimensionless total free energies TFE 
versus circle-circle separations «r for the case where inverse 
screening length times circle radius Ka—1.8, and as a function 
of circle edge-line separation rec. Note that in order to obtain a 
three dimensional free energy one must multiply TFE by K g | T 
and by the number of Bjerrum lengths the circular cylinders 
extend into the z direction. The circle boundary condition 
(b.c.) is .0507 constant potential slope (corresponding to a -1 
free surface potential if «a=1.8) and the confining line b.c. is 
-2 constant potential. Lines are two parameter fits of the data 
points to Eq. 2, thereby obtaining the geometric factor GF 
and 

^local *^ 

function of circle edge-line separation kc. The 
repulsive interaction is substantially reduced and the local 
Debye screening length is shortened by the presence of the 
charged plates. 

Figure 2 shows seven EDL potential barriers between 



two circles confined by line charges arranged in the ge- 
ometry depicted in Fig. 4. The dimensionless total free 
energy TFE (Eqs.5-9) is plotted versus the circle's cen- 
ter to center separation r. The circles have a free po- 
tential equal to -1, but it is their surface charges which 
are held constant during the circle-circle interaction, this 
because the circle surfaces are assumed to posess strong 
acid groups which resist charge condensation. Their con- 
stant potential slopes are equal to .0507, the appropriate 
slope for a -1 potential free surface when «a=1.8, i.e. the 
screening length is 20 grids, and the circle's radius a is 
fixed at 36 grids. Here the confining line charges are held 
at constant -2. potential, and are thus said to be perfectly 
regulating because their surfaces are assumed to be com- 
posed out of weak acid groups. Circle edge-line charge 
separation c is varied from \ to 6 kT 1 . Circle motion 
is assumed to be adiabatically cut off from the motion 
of the ion gas, i.e. the ions are assumed to readjust to 
the new circle configurations with extreme rapidity. Cir- 
cle separation r begins at 15k -1 and is decremented by 
units of .5 k _1 down to 9 re -1 , all the while solving the 
PB equation and subsequently recording the configura- 
tion's Hclmholtz free energy TFE. A two parameter fit 
of the (r,TFE) data points is made to the 2d EDL po- 
tential function (Eq.2), thereby obtaining the interaction 
screening length i^ ooal , and observing that this circle- 
circle interaction screening length is decreased from the 
bulk value by close confinement of the line charges. 

Figure 3 breaks down the relative contributions made 
to the circle-circle interaction potential between the elec- 
trical, entropic, and chemical free energies (Eqs.5-9) for 
Fig.2's case of strongly confined circles. The case of two 
constant potential line charges confining two constant 
charge circles contains two attractive and one repulsive 
terms. The smaller of the attractive terms is the chemical 
free energy term (the number of line-charge counterions 
increasing as the circles move together). The larger at- 
tractive term is electrical in origin arising from the coun- 
terions located in between the two circles — like in a hy- 
drogen molecule. However, at 300 Kelvin, the repulsive 
entropic term dominates the attractive terms. 
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FIG. 3. The three component terms of the circle-circle in- 
teraction free energy for Fig.2's case of strong confinement are 
plotted versus circle center-center separation Kr, where EE is 
electrical energy, always attractive, — TAS due to ion/solvent 
configurational entropy, always repulsive, and CFE is a weak 
attractive chemical free energy (line charge counterion num- 
ber increases when circles are close together). Notice that the 
repulsive ion/solvent configurational entropy term is domi- 
nant, meaning that the total potential is repulsive and that 
this repulsion is due to osmotic pressure forces. The condi- 
tions are two constant -2 potential line charges strongly con- 
fining (kc=.5) two interacting constant potential slope circles 
(.0507 slope appropriate to a free -1 potential if ko, = 1.8). 

Figure 4 shows three contour plots associated with 
Fig. 2 and appropriate to greater or lesser degrees of con- 
finement. When the circles are more confined the con- 
tours surrounding each circle become more elongated, re- 
flecting the fact that many of the circle's electric field 
lines are now terminating locally in the line charge's dou- 
ble layer. One could think as if the field lines were termi- 
nating in partial image charges created in the line charge 
double layer. 




FIG. 4. Contour plots of the dimensionless potential (<f>) 
for the conditions shown in Fig. 2. The circles have constant 
charge (their boundary potential slope which equals .0507 cor- 
responds to an unconfined potential equal to -1 when Ka— 1.8) 
and the line charges are at constant -2 potential. All poten- 
tials are negative. jcr=10 for all three plots. There is no 
electric field within the black circle interiors, (a) The cir- 
cle-circle interaction is virtually undisturbed by the distant 
charged lines, kc~6. (b) The circle-circle interaction is be- 
coming weakened by confinement, kc=2. Note the increasing 
circle surface potential, (c) Strong confinement, kc=.5. Field 
lines normally involved in the circle-circle repulsive interac- 
tion are being redirected, terminating on the counterions in 
the charged line's double layer. 

Figures 5,6, and 7 are all connected. In Fig. 5 the con- 
fining line potential is increased to -6 and the confined 
circle-circle interaction screening length Kj~ ocal is appre- 
ciably lessened when compared with Fig. 2 which has only 
a -2 potential confining line. The large line capacitance is 
acting to hide the circles from each other. Figure 6 gives 
the relationship between confining line potential and lo- 
cal screening length, while Fig. 7 shows for a free line 
charge the relationship between confining line potential 
4>iine and line capacitance per unit length C/L. 

C/L = V , d " ' (10) 

01ino 
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FIG. 5. Reduction of interparticle repulsion by increas- 
ing the confining potentials compared to Fig. 2. Plots of total 
free energy TFE versus circle-circle separation (/tr) for .0507 
constant potential slope circles confined by line charges held 
at a -6 potential, and k<x=1.8. The circle-circle barrier po- 
tential is reduced by approximately one order of magnitude 
when the circles are severely confined. Note that the degree 
of confinement nc is varied from 6 to i . 



like inverse screening length n for a free line when 4>une is 
held at zero potential (obtained by dividing the r.h.s. of 
the Grahame equation by 4>u ne , while taking the limit as 
0iinc — > 0). In Fig. 6 we plotted local screening length ver- 
sus confining line potential and obtained gaussian func- 
tions symmetric about the zero of potential. Notice that 
the local screening length is .85k -1 for zero potential 
line charges strongly confining the circle-circle interac- 
tion, kc=.5, and that this is connected with the fact 
that a zero potential confining line charge posesses an in- 
creased capacitance determined by a Grahame equation 
modified to account for the confinement jl6| . Essentially, 
strong confinement brings the zero potential line an ad- 
ditional capacitance due to its sharing of the confined 
circle's counterions. Interestingly, the local interaction 
screening length is independent of the sign of the surface 
charge on the confining lines. Apparently as regards the 
local screening length it does not matter if a negative 
field line sourced from a circle terminates in the line's 
double layer on a positive counterion or on the positive 
line charge itself. 

1.5 i ■ 1 ■ 1 ■ 1 ■ 1 
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FIG. 6. The ratio of local screening length K ; ~* ai to bulk 
screening length k _1 as a function of the confining line poten- 
tial 4>nne boundary condition. The curve is symmetric with 
respect to the confining line potential. The conditions are 
na—1.8, solid circles correspond to nc = .5 (strong confine- 
ment), hollow squares to kc = 1 (weaker confinement), while 
the circles have constant (charge) potential slope =.0507. 

Figure 7's potential slope fits the Grahame equation 

^ = K 2sinb^. (11) 
dn 2 v 1 

Also note Fig. 7 capacitance per unit length's nonzero y- 

intercept (~ .05 = n) shows that line capacitance goes 
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FIG. 7. Plots for a free line charge of its boundary 
potential slope <r(= dl ^H° c , solid circles) and its intrinsic ca- 
pacitance per unit length C/L (= 17 , hollow squares; see 

Eq.10) versus free line potential cj>u ne . The line tracing the 
solid circles is a fit to the Grahame equation (Eq.ll). 

The second Eq.2 fitting parameter GF stands for the 
DLVO geometric factor |f| , which accounts for the inter- 
nal volume of the cylinder being excluded to the screen- 
ing ions. For a typical colloidal particle there are two 
effects competing to set effective surface charge Z\J GF 
as one increases the concentration of potential determin- 
ing ions (p.d.i.). A p.d.i. electrolyte is contrasted to an 
indifferent electrolyte, the distinction being that an in- 
different electrolyte is unable to bind to surface residues. 
The first effect concerns binding of a p.d.i. to a surface 
residue, thereby decreasing Z, while the second occurs 
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as the concentration gradient of p.d.i.'s in the vicinity 
of a charged surface is increased, thus making the par- 
ticle appear to have increased its surface charge. This 
second effect is the one accounted for by DLVO's GF. 
If one holds Z constant while increasing ionic strength, 
then one observes an increasing effective particle surface 
charge Z* = Z\J GF due to the steepened counterion con- 
centration gradient. We see in Fig. 8 that GF increases 
in a gaussian fashion with increasingly negative confining 
line potential. Compare the behaviors of the two fitting 
parameters GF (Fig. 8) and n~[^ cal (Fig. 6). Local screen- 
ing length doesn't care about the sign of the surrounding 
line charge, but GF certainly does. 
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FIG. 8. Plot of GF = versus line potential 
4>une for two confined interacting .0507 constant potential 
slope circles when K,a—l.S. is the effective surface charge 
on the circle. The solid circular data points correspond to 
kc = .5 (strong confinement), and the hollow squares to 
kc = 1 (weaker confinement). 

After K;~^a; as a function of c is obtained from Eq.2 
(see Figs. 2 and 5) we make a one parameter fit to the 
function 



FIG. 9. Local screening length versus the confinement 
distance between the circle and the line charge, kc, for various 
circle and line boundary conditions, where <fr = K ^ T / e is a di- 
mensionless surface potential boundary condition (b.c.) and 
a = |£ a surface charge b.c, and with Ka = 1.8. The lines are 
one parameter fits to the factor A in Eq.12. The solid data 
points demonstrate the line potential's control over the local 
screening length. Meanwhile changing the circle b.c. had little 
effect on interaction screening length, as seen by comparing 
the solid circle, open circle, and open square data sets. Note 
that a constant charge b.c. on the confining plate correspond- 
ing to a -2 potential, i.e. the open diamond (^=.270) had a 
significantly smaller effect on local screening length than did 
its -2 confining potential counterpart. This probably had to 
do with the fact that line charge counterion number increased 
when the circles were close together, for the case of constant 
potential confining lines (see Fig.3). 

Fig. 10 shows that decreasing c, i.e. increasing the 
degree of confinement, increases the GF for all of the 
constant-charge circle-circle interactions. In contrast, for 
the case of confined constant-potential circles the square 
data points show that the combination jjGF actually 
decreases. For strong confinement (small c) condensa- 
tion of counterions decreases the circle's surface charge 
^, and this effect wins out over an increasing GF. 



local 



(12) 



and plot the result in Fig. 9. Local screening length 
goes like the inverse square root of local ionic strength 
( K k>cal ~ n ' 5 )j which itself goes like an exponential func- 
tion of the local potential (Boltzmann factor: n <~ e~^) 1 
and finally this local potential goes like a decaying ex- 
ponential function of distance from a charged surface 
(0- e- KC ). 



6 



CD 



N 



20000 



15000 



10000 - 



5000 



circle b.c. 


line b.c. 


\ 




c $ =-1 


=-2 


□ a =0.0507 


a=0.1131((|)=-2) 


\V ■ a =0.0507 

\^^a^a05O^ 


<|> =-4 
<|) =-6 


_a « Q — 







2.5 



7.5 



KC 



FIG. 10. Squared surface charge per unit length times 

the DLVO geometric factor GF (see Eq.2) versus con- 
finement kc for various circle and line boundary conditions, 
where <j> — K] ^ T / e is a dimensionless surface potential b.c, 

a — 4$ a surface charge b.c, and na — 1.8. Only in the 
constant potential circle case where surface charge regulation 
acts to decrease Z does GF decrease with increasing con- 
finement. 

Finally we vary bulk screening length from 6 grids to 30 
grids for various amounts of confinement, and obtain the 
dependence of local screening length on bulk screening 
length for these differing degrees of confinement (Fig. 11). 
In our 2d geometry there are 3 independent length scales : 
circle radius (fixed at 36 grids), circle edge- line separation 
c (varying from 10 to 40 grids), and k _1 (varying from 
6 to 30 grids, i.e. na decreases from 6 to 1.2). The 
circle-circle interaction's length scale Unocal 1S a function 
of these three. As one increases the bulk screening length 
one decreases the space's midplane potential towards the 
line potential (cf)ii ne — —2 in Fig. 11), and thus more of the 
circle's field lines traverse a region where the local ionic 
strength has been increased by a factor of cosh[— 2] = 
3.76 times the bulk ionic strength (counterions increased 
by the factor e 2 while the coions decrease by e~ 2 ). Hence 
the screening length "seen" by these field lines will be 
reduced by a factor ^ = .52. While for the case k -1 = 
6 most of the interaction field lines sample electrolyte 
where k 
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FIG. 11. Local screening length function of 

bulk screening length for various degrees of confinement. The 
circles have a constant potential slope =.0507 (corresponding 
to a -1 free surface potential only when na = 1.8) and the line 
charges are held at a -2 potential. Local screening lengths 
were obtained by a fit to Eq.2 of the total free energies in a 
region where circle surface to surface separation was approxi- 
mately 5.5 bulk screening lengths. The diagonal corresponds 
to ^Yocal = °P en circles to c=10 grids, open squares to 

c=20 grids, open diamonds to c=40 grids, while the circle 
radius is 36 grids. The local screening length increases to 
the bulk screening length as the confining line charges move 
apart. 

Solutions to the confined non-linear Poisson- 
Boltzmann equation for like- charged cylinders appear 
to be repulsive, screened, Yukawa-like potentials, having 
an effective screening length n^ cal which is found to be 
a decreasing function of increasing confinement. The ef- 
fective charge of the cylinders Z* = j-^/GF is found to 
be an increasing function of increasing confinement. 

III. CONCLUSIONS 

By numerically solving the two dimensional, non- 
linear, confined Poisson-Boltmann equation we've found 
that the electrical double layer (EDL) repulsion can be 
significantly decreased between two parallel cylinders 
when confined between two charged plates because the 
confining double layers constitute a high capacitance re- 
gion which screens the cylinder-cylinder interaction. 
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